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Abstract: A technique based on using optical fiber taper waveguides for 
probing single emitters embedded in thin dielectric membranes is assessed 
through numerical simulations. For an appropriate membrane geometry, 
photoluminescence collection efficiencies in excess of 10 % are predicted, 
exceeding the efficiency of standard free-space collection by an order of 
magnitude. Our results indicate that these fiber taper waveguides offer 
excellent prospects for performing efficient spectroscopy of single emitters 
embedded in thin films, such as a single self-assembled quantum dot in a 
semiconductor membrane. 
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1. Introduction 

The development of novel techniques for efficient detection and spectroscopy of individual 
quantum emitters in the solid state is essential for an understanding of the emitter and its re- 
lationship with the surrounding environment. Challenges involved with spectroscopy of such 
systems include tight focusing requirements and separation of the signal of interest from back- 
ground scattered light due to the host crystal. In recent years, tools such as high numerical 
aperture objectives |[Tl|2][3l and near-field scanning optical microscopy tips ||4| have been used 
for both efficient photoluminescence (PL) collection and resonant optical spectroscopy. 

A related problem is optical spectroscopy of microphotonic resonators, where wavelength - 
scale focusing and effective separation of the signal from background scattered light are also 
necessary. The optical fiber taper waveguide is a tool that has been investigated for such exper- 
iments. These structures, sometimes referred to as silica nanofibers, are standard single mode 
optical fibers that have been heated and stretched down to a wavelength-scale minimum diam- 
eter At such small dimension, the evanescent field of the waveguide mode extends into the air 



cladding and can be used to interrogate surrounding structures. By conducting the tapering pro- 
cess adiabatically and symmetrically, a double-ended device with standard single mode fiber 
input and output can be produced, with a typical overall transmission loss of less than 10 %. 
The combination of low loss, single mode guidance, and access to a wavelength-scale evanes- 
cent field has made tapered fiber probes invaluable for a wide number of microcavity-based 
experiments, including those in single atom and quantum dot cavity quantum electrodynam- 
ics ||5]|6]|7l. These properties also suggest that they can be effective tools for single emitter 
spectroscopy. Indeed, theoretical work mm has indicated that silica nanofibers can be used to 
efficiently collect fluorescence from single gas phase atoms, and experimental progress to this 
end has been made ifTOl . 

The above studies indicate that optical fiber taper waveguides can be a general spectroscopic 
tool for nanophotonic systems, providing motivation for the present work. Here, we study their 
potential for efficient non-resonant and resonant fluorescence collection experiments on single 
solid state emitters embedded in a dielectric slab. Additional motivation is provided in previous 
experimental work fTTl, where a fiber taper waveguide was used in non-resonant PL measure- 
ments of a single InAs quantum dot embedded in a 256 nm thick GaAs membrane on top of 
an Alo.7Gao.3As pedestal. There, a luminescence collection efficiency of w 0.1 % was esti- 
mated based on the measured saturated photon count rates, neglecting possible radiation rate 
modification due to the host structure, a non-resonant microdisk cavity, and any quantum dot 
non-idealities. Important questions left unanswered include how close this efficiency was to the 
theoretical maximum for the particular geometry considered, and how the geometry may be 
modified to lead to higher efficiencies. Furthermore, the possibility of studying resonant fluo- 
rescence, not addressed in ifTTl . should also be considered. As described in other recent works 
m El El El El lT2ll - resonant spectroscopy has advantages in comparison to non-resonant PL 
measurements in terms of greatly improved spectral resolution, potentially improved temporal 
resolution, the ability to avoid generation of decoherence-inducing excess carriers, well-defined 
state preparation, and possible utility in quantum information processing experiments. 

To address these points, we consider the spontaneous emission of a dipole embedded in a 
semiconductor slab surrounded by air and in contact with a micron diameter silica optical fiber 
Specifically, we envision an experimental measurement setup as depicted in Fig. El where the 
optical fiber taper waveguide is brought into close proximity with a thin semiconductor mem- 
brane hosting a single quantum emitter The fiber has been adiabatically tapered to reach the 
single-mode condition over its central region (hundreds of microns long). The emitter may be 
optically pumped through an external fiber-coupled source, and radiates in both directions into 
the fiber, with a total coupling efficiency T7pl into the fundamental fiber mode (both backward 
and forward directions). We perform full wave, 3D finite difference time domain (FDTD) elec- 
tromagnetic simulations of a single classical dipole radiating in the probing structure shown in 
Figs [nb)Etc), and obtain the fiber-collected power We find that fiber collection efficiencies in 
excess of « 12 %(30 %) may be achieved for dipole moments horizontally (vertically) oriented 
with respect to the membrane plane, exceeding what may be achieved with standard free-space 
optics PL collection schemes by an order of magnitude or more. These results suggest that op- 
tical fiber taper waveguides are potentially quite valuable in future studies of single solid state 
quantum emitters. 

To determine ways of further improving collection efficiency, we analyze the dependence 
of both the modified spontaneous emission rate and PL collection efficiency on the electro- 
magnetic supermodes of the probing structure. The starting point of our analysis is Jsl, which 
determines the modified spontaneous emission rate of a multi-level atom near a silica fiber 
Addition of the semiconductor slab is a non-trivial modification of this problem, requiring the 
use of finite-element method (FEM) simulations to compute the guided and radiation modes 




Fig. 1 . (a) Single emitter probing setup based on a tapered fiber waveguide, (b) Detail of (a), 
showing the membrane that carries the emitter, (c) Schematic of the substrate cross-section, 
showing membrane and fiber. 



of the composite slab-fiber system. These modes are used along with FDTD to determine the 
individual supermode contributions to the spontaneous emission fiber coupling fraction T7pl, 
under the assumption that the embedded emitter can be approximated as a two-level system. 

The paper is organized as follows. In Section |2] we use FDTD to estimate the spontaneous 
emission collection efficiency into an optical fiber for a single quantum dot embedded in a 
semiconductor membrane, as shown in Fig.[T] The results of this section are analyzed in Sec- 
tion[3]in terms of the propagating supermodes of the probing structure. In SectionH] we further 
discuss the physical interpretation of these results and consider how they change with fiber size 
and if the host semiconductor membrane rests on a substrate (i.e., a non-undercut slab geome- 
try). We also consider a pair of specific configurations with improved collection efficiencies, for 
which additional investigations are underway. Finally, we describe how resonant fluorescence 
measurements can be performed through our fiber-based probing and collection scheme. 

2. Fiber-based embedded single emitter photoluminescence collection 

In this section, we analyze the photoluminescence (PL) collection from a single emitter embed- 
ded in a dielectric slab with the fiber taper probe setup illustrated in Fig. [T] For non-resonant 
PL measurements, application of the fiber probe in the configuration of Fig.|2ja) is envisioned, 
which provides fiber-confined paths for excitation and PL collection in both forward and back- 
ward directions. As indicated in Fig.[Tlb), the micron-scale single mode region of the optical 
fiber waveguide is brought into close proximity with the top surface of a dielectric membrane 
that hosts the emitter, over a length of several to hundreds of wavelengths. Fiber and membrane 
together form the composite dielectric waveguide with cross-section shown in Fig. [lie). The 
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Fig. 2. (a) Envisioned experimental configuration for fiber-based non-resonant photolumi- 
nescence (PL) spectroscopy with a tapered fiber waveguide. Emitted fight is coupled into 
both forward and backward channels of the fiber taper waveguide, and can be wavelength 
resolved with a grating spectrometer or spectral filter, (b) Schematic of single emitter exci- 
tation and PL collection via the tapered fiber probe. A non-resonant pump signal is injected 
into the input fiber and converted into a guided supermode of the composite waveguide, 
illuminating the slab-embedded dipole. The dipole radiates into guided and radiative su- 
permodes, with rates Fq and Fy, respectively. Power is transferred with efficiency /o from 
the supermode to the fiber mode and vice-versa. 



composite waveguide supports a complete set of guided, leaky and radiation supermodes ifTSl . 
originating from the hybridization of fiber and slab modes. As depicted in Fig. |2lb), part of 
the non-resonant pump power initially carried by the fiber is coupled to a guided supermode of 
the composite waveguide and transferred toward the emitter Illuminated by the pump light, the 
emitter radiates (at a red-shifted wavelength) into modes of the composite fiber-slab waveguide. 
Our objective is to determine the fraction of this emission that would be measured at the single 
mode fiber output, after the micron-scale fiber has transitioned away from the slab. The portion 
of the total emission collected by the fiber is carried by those supermodes with sufficient trans- 
verse confinement to reach the fiber transition regions with significant amplitude. Thus, apart 
from guided modes, leaky modes may contribute significantly to the total PL collection. 

Within this scenario, we estimate the spontaneous emission collection efficiency for a sin- 
gle emitter in a geometry that approximates the tapered-fiber-based quantum dot spectroscopy 
setup of IITTII . In that work, PL spectroscopy of single self-assembled InAs quantum dots em- 
bedded in a non-resonant 256 nm thick GaAs microdisk was demonstrated, with an estimated 
collection efficiency of 0.1 % (including both channels of the optical fiber). Our numerical 
estimate here is obtained through full electromagnetic simulations of a point dipole radiating 
inside the composite waveguide of Fig.[TJc). The geometry consists of a dielectric slab of index 
"slab = 3.406, and fiber of radius r = 500 nm and refractive index «fiber = 1-45 placed directly 
on the top surface of the slab. The dipole is assumed to radiate at A = 1.3 \lva, a wavelength at 
which the fiber possesses, apart from the fundamental guided mode, a near-cutoff mode which 
is not significant in our calculations. Slab membranes of thicknesses varying between 100 nm 



and 260 nm were investigated, sufficient to avoid quantum dot dephasing due to surface-state- 
related processes, as suggested in |[T4l. In SectionHl we consider how these resuhs change when 
r is decreased to 300 nm and when the slab membrane is placed directly on a substrate. 

2.1. Simulation Method 

To estimate the PL collection efficiency of our fiber-based probing scheme, we considered a 
single classical electric dipole radiating in the composite dielectric waveguide of Fig. [He). An 
upper bound for the percentage of the total emitted power Prot. coupled to the fundamental 
optical fiber mode at an arbitrary position z along the guide is 

p 

J7PL = 2— ^/fiber, (1) 
^ot. 

where is the power flowing normally through the constant-z plane, and /fiber is the overlap 
integral IHl 

Re {x/,(e/ X h*) • zdS //^(e x h)) • zds] 

/fiber — T ^ • (2) 

Re [jJsitf X h}) ■ zdS\ Re {//.(e x h*) • zdS} 

Here, {e,h} are the steady-state radiated fields at position z and {e/,h/^} the fundamental fiber 
mode fields. The factor of 2 in Eq. ([l) is introduced to account for collection from both fiber 
ends. 

The steady-state fields were obtained through FDTD simulations lfT6l . The symmetry of 
the geometry allowed us to choose either anti-symmetric (x x E 0) or symmetric (x x H = 
0) boundary conditions on the yz-plane, the former being used for x-polarized dipoles, and 
the latter for y- and z-polarized dipoles. The computational domain was cubic with a 10 /im 
side and perfectly-matched layers (PMLs) were used around the domain limits to simulate an 
open domain. Dipole excitation consisted of a X = 1.3/im carrier modulated with a 0.1 nm 
bandwidth gaussian envelope. Simulations ran until no field amplitude could be detected in 
the domain. The total radiated power Prot. was calculated by adding the integrated steady- 
state power through all the computational window sides. The PL collection efficiency tjpl was 
calculated at several z-planes along the propagation direction (z) in the computational domain. 

2.2. Simulation Results 
2.2. 1 . Total radiated power 

Figure Oa) shows the FDTD-calculated, total radiated power (Prot.) for x-, y- and z-oriented 
dipoles located at the slab center in the structure of Fig. [Tfc). The curves were normalized 
to the total radiated power of an electric dipole current source J(r,f) = — ia)p5(r)e^''"o' in a 
homogeneous dielectric medium of refractive index «siab and electric dipole moment p ifTTl : 

^om. = ^«slab|pp|^. (3) 

The total radiated power is not significantly different when the fiber is absent. This is ev- 
ident in Fig. Ob)- where the power radiated by vertical or horizontal dipoles into modes of 
an isolated slab are plotted. These curves were calculated using the transfer-matrix method 
of ifTSl . Transverse-electric waves (TE) have electric field components in the slab plane only, 
while transverse-magnetic (TM) modes have only in-plane magnetic field components. As a 
result, a vertical (y-directed) dipole excites only TM waves, while horizontal {x- and z-directed) 
dipoles can excite both TE and TM waves lITSl . although as we see in Fig. Ob), the TE mode 
excitation is much more significant. The similarity between the total radiated power curves 
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Fig. 3. (a) Total radiated power at X = 1.3/im for a dipole in the composite slab-fiber 
waveguide, normalized to the radiated power in a homogeneous medium of index «siab- 
(b) Guided and total radiated powers into TE (TM) waves for a horizontal (vertical and 
horizontal) dipole at the center of a dielectric slab with /ijiab = 3.406. 



for the slab with or without the fiber indicates its weak perturbative effect on the emission, 
so we first proceed by analyzing the isolated slab. Figure |4{a) displays the effective indices 
of the TEo, TEi and TMq modes, for varying slab thickness. The evolution of the field con- 
centration at the dipole for both TE and TM modes may be inferred from the effective length, 
Lgff = / e|Ept/y/|E(ro)|, where ro is the dipole position, plotted in Fig. Hfb). The effective 
length is inversely proportional to the dipole power coupled to the mode lfT9l |20| . 

Given the large refractive index of the slab relative to air, total internal reflection leads to 
preferential dipole excitation of guided waves. This is verified in the curves for guided and total 
radiated powers in Fig.[3jb) il8] . For a horizontal dipole at the slab center, most of the radiated 
power is carried by the TEq mode, which has a relatively small effective length, Leff- The TEi 
mode is not excited, since its electric field has a node at the dipole location. The TMq mode is 
also not excited, as E, = and E- has a null at the slab center. For a vertical dipole, excitation of 
the TMq mode is more effective, especially since no TE modes are generated ifTSll . The strong 
emission inhibition for vertical dipoles at small t^]ah is due to poor outcoupling into the air, 
resulting from total internal reflection, and poor field concentration, apparent in the large TMq 
effective lengths Leff of Fig.|4tb). As a result, even if high fiber collection efficiencies can be 
achieved, this strong suppression of spontaneous emission can lead to low overall photon count 
rates at the detector. 

2.2.2. Photoluminescence Collection 

Using Eq. ([l]), the collection efficiency T7pl was calculated for a range of propagation lengths 
a few wavelengths away from the dipole location (1 fim < z < 5 fim). The results reveal strong 
periodic oscillations as a function of this distance, as shown in the inset of Fig. |3 a). As we shall 
discuss in the following section, this oscillation is due to the beating of two principal waveg- 
uide supermodes with distinct propagation constants, corresponding to the sloshing of power 
between fiber and slab, as commonly observed in coupled waveguide systems fT5\. Oscilla- 
tions in TjpL are accompanied by a decrease in its peak value - as we shall see in Section |3] this 
decay is due to leaky modes with propagations lengths of several wavelengths. The maximum 
and minimum values of T7pl within the simulated propagation range were taken to generate 
the PL efficiency plots for x-, y- and z-directed dipoles of Fig. |5] Note that the values plotted 
in Fig. |5] include collection through both the forward and backward channels of the fiber, an 
experimentally viable situation through use of a fused fiber coupler (Fig.|2a)). 
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Fig. 4. (a) Effective indices of tiie slab TEq, TEi and TMq modes at A = 1.3 /im. The HEn 
fiber mode effective index is also shown. Inset: field components of TEq, TEj and TMq 
slab modes, (b) Effective lengths, Lgff, for the TEq and TMq slab modes. 



We first consider the results for in-plane oriented dipoles (in the x and z directions), which are 
of relevance to quantum dot experiments such as in ifTDl . For jc-oriented dipoles, the maximum 
collection efficiency ranges from « 7 % to w 2 % for t^xab varying from lOOnm to 260 nm, 
while the minimum efficiencies remain close to 0. 1 % throughout. For z-oriented dipoles, the 
efficiency varies below 1 % for all thicknesses. In ifTTI . a single quantum dot inside a 256- nm 
GaAs slab was probed, and a measured collection efficiency of 0. 1 % was estimated. This value 
is within the limits predicted by Figs.|5ja) and|5jc), given that the in-plane dipole orientation 
was not known, and the quantum dot was assumed to have unity radiative efficiency. In addition, 
the in-plane position of the quantum dot was such that the underlying Alo.7Gao.3As support 
pedestal could have had significant influence on the PL collection. We consider its effects in 
Sectiongl 

To further gauge our current results, we use the method of ifTSl to plot, in Fig.|5j the max- 
imum achievable free-space collection efficiency for a dipole embedded in the slab (without 
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Fig. 5. Fiber-collected PL coupling efficiency 7]pl for (a) x-, (b) y- and (c) z- oriented 
dipoles at A = 1.3 /xm. Collection through both the forward and backward fiber chan- 
nels is considered (notice changing axis scales). Since T]pl oscillates with z, its maximum 
and minimum along z <5 /im are plotted. The inset in (a) shows T7pl as function of z. 
Squares: FDTD results. Circles: Mode-expansion results. Solid lines: Free-space collection 
efficiency considering a NA=0.7 objective. 



the fiber), using a 0.7 numerical aperture (NA) objective. We see that the maximum achievable 
T7pl for horizontally-oriented dipoles in the fiber-based collection scheme is superior to the 
free-space collection efficiency by as much as an order of magnitude, as indicated in Fig.|5ja). 
As is the case with other co-propagating waveguide directional couplers ifTSl . control of the 
interaction length is necessary to achieve the most efficient power transfer. Here, to obtain the 
highest values of tjpl, one might envision controlling the length of the membrane through stan- 
dard planar fabrication techniques. Alternately, the effective fiber-membrane interaction region 
can be limited using 'dimpled' fiber taper waveguides II2TI . Given the low propagation losses 
in the optical fiber and a typical insertion loss of 10 % through the taper, the present technique 
constitutes a competitive option for PL spectroscopy. 

For vertically oriented dipoles, the collection efficiency may reach above 35 % at fjiab = 
140 nm, where the dipole emission is strongly suppressed. For thicker slabs, the suppression is 
reduced, but the collection efficiency drops, reaching « 3 % for fsiab = 250 nm. Even this col- 
lection efficiency is nonetheless an order of magnitude higher than what may be obtained with 
free-space collection with a 0.7-NA objective, plotted in Fig. Elb). Although self-assembled 
quantum dots in semiconductor slabs are typically in-plane oriented dipoles, we note that ma- 
nipulation of the quantum dot growth can change its shape sufficiently to produce structures 
that support TM-polarized emission 1221 . The results presented here would be of significant 
interest to such systems. 

While the results presented in this section are analyzed in detail in Section[3] we provide here 
a brief explanation for the relatively high collection efficiencies found, regarding the fiber-slab 
system as a directional coupler. In a directional coupler, power transfer takes place between 
waves, in respective waveguides, of similar propagation constants. Although transfer is most 
efficient when both waves are perfectly matched, small but significant transfer may still oc- 
cur between mismatched modes. As discussed in Section l2.2.1l most of the x-polarized dipole 
radiation is carried by TE slab mode waves. A significant portion of these waves, though phase- 
mismatched to the fiber mode, transfer small amounts of power to the latter, collectively leading 
to relatively high (8 % maximum) collection efficiencies. The substantially higher (40 % maxi- 
mum) collection efficiency obtained for y-oriented dipoles results from a more effective power 
transfer from TM slab waves to the fiber mode: phase matching between the fiber and TMq slab 
modes is apparent in Fig.Ua), for slab thicknesses close to 160nm. 

3. Analysis 

While the FDTD simulations provide estimates of the total spontaneous emission rate and fiber- 
coupled collection efficiency, they do not necessarily provide physical understanding of effects 
such as the oscillation and decay in PL collection efficiency as a function of separation from 
the dipole in the propagation direction, or the insight needed to extend this work to modified 
geometries. In this section, we address these issues through FEM simulations that determine the 
modal structure of the composite fiber-slab system. This allows us to determine the spontaneous 
emission rate, and contribution to the total collected PL, of individual supermodes of the sys- 
tem. This approach furthermore allows us to estimate the evolution of the coupling efficiency 
for coupling lengths much longer than practical in FDTD simulations. 

3.1. Supermodes of the composite waveguide 

The composite waveguide of Fig. [T] supports a set of hybrid supermodes of the fundamental, 
micron-scale fiber mode and the bound, TE,„ (£,. = 0) or TM,„ {Hy = 0) slab modes. These su- 
permodes are henceforth referred to as hybrid-TE,„ or hybrid-TM,„. The fiber creates a region 
of increased effective refractive index on the slab plane, spanning a discrete set of laterally con- 
fined hybrid-TE or hybrid-TM supermodes. A continuum of radiation hybrid-TE or hybrid-TM 



supermodes also exists, which represents waves not laterally confined by the increased fiber- 
induced effective index. A portion of this radiative supermode spectrum corresponds to optical 
power that, although unconfined, lingers near the fiber for relatively long propagation distances 
(tens of wavelengths), and may contribute significantly to the PL collection efficiency. These 
wave bundles may be represented by a discrete set of leaky supermodes, which extend infinitely 
in the lateral direction and decay in the propagation direction (TSl- The relevant waves for PL 
collection within our scheme are thus the guided and leaky hybrid-TE or hybrid-TM super- 
modes. 

The supermode field distributions of our structure were obtained through the vector eigen- 
value problem for the magnetic field H 



Here, H = H(r) exp(^z), and £, = — = icoiies/c is a complex propagation constant with a 
phase term )3 and a decay term a,, or equivalently a complex effective index ngff- The number 
of free-space wavelengths necessary for the supermode amplitude to decay by a factor < 5 < 
I is Ns = — n,^' ln(5)/2;r, where n^^ = Im{/ieff} This number is used in our analysis to 
roughly determine the probing length (i.e., the length over which the fiber contacts the slab) for 
supermode power collection to be most effective. Since the dipole is aligned with the center of 
the fiber, the problem is symmetric with respect to the yz plane. An x-oriented dipole thus only 
excites hybrid-TE supermodes, for which x x E = on the symmetry plane. Conversely, y- and 
z-oriented dipoles only excite hybrid-TM supermodes, for which x x H = on the yz plane. 

The isolated slab waveguide supports a TEq and a TMq mode for all slab thicknesses fsiab 
considered, and a TEi mode for fsiab ~ 200nm (Fig.lHi, so that the composite waveguide struc- 
fiire supports hybrid-TEg and hybrid-TMo and hybrid-TEj modes in the same ranges. Field 
profiles for the strongest-confinement hybrid-TEg and hybrid-TMo supermodes are shown in 
Figs. Ha) andlStd) of Section[331 

3.2. Modified spontaneous emission rate in the composite waveguide 

To describe the spontaneous emission of a single quantum emitter into the composite dielectric 
waveguide of Fig. [T]in terms of its supermodes, we begin by making use of IS), where the 
spontaneous emission rate of a multi-level atom into the guided and radiative modes of a silica 
optical fiber are calculated. These expressions are derived from the Heisenberg equations for 
a multilevel atom experiencing a dipole-type interaction with a vacuum field reservoir that is 
described in terms of propagating waveguide modes. While [ 8 | focused exclusively on optical 
fibers, here we deal with general dielectric waveguides, so that the expressions we use are 
modified accordingly from their originals, with the derivation presented in Appendix B. 

The waveguide, extending in the z direction, is open, and thus supports a finite set of guided 
modes and a continuum of unconfined, radiation modes lfT3l . Guided modes and radiation 
modes are respectively labeled with the indices /i and v, and have propagation constants 
and j3 at frequency o. For the radiation mode continuum, j3 may assume any value in the in- 
terval |j3| < (On/c, where « is the refractive index of the medium surrounding the waveguide. 
Evanescent radiation modes are excluded, as they do not participate in the radiation process 



For a two-level atom located in a general dielectric waveguide at a position ro, with a tran- 
sition of energy Tkoq and dipole moment p, the spontaneous emission rates into guided and 
radiative modes are, respectively, = 27r|Gn(a)o)P and Fy = 2%\Gv{(0q)\'^ , where 




(4) 



Here, and ey are the electric fields of the guided and radiative modes, with respective prop- 
agation constants and j3. Also, = d^^^/da and A^y is such that |fT3l 

JJdS{tyXhl,)-z = Nydij5-l5'). (6) 
s 

Guided and radiative modes are normalized so that 

||^/5e(r)|e^|2 = l (7) 

s 

jj dSe{r){ey ■el,)p=p>,„=„, = d{co- co') (8) 
s 

where S is the cross-section of the waveguide and e(r) the spatial permittivity distribution. The 
total emission rate is 

r = Lr^+E/^i8rv. (9) 

H n •' 

According to Eq. dA-Sb . the power emitted from a classical dipole into a particular mode 
depends on the modal electric field concentration at the dipole position (the modal electric 
fields are normalized to the power density flux integrated over the waveguide cross-section). We 
introduce the effective area Aeff = / e|Epc/5/|E(ro)P as a measure of the field concentration 
at the emitter position. This expression is implicit in the modal rates of Eq. (|5]l, given the 
normalization expressions, Eqs. d?) and (O. The total spontaneous emission rate, being the 
sum of the individual modal rates, deviates considerably from the free-space rate if highly 
concentrated modes (i.e., small effective areas) exist at the emitter location. 

To avoid calculation of all radiative modes of the composite waveguide, the spontaneous 
emission rate Y was obtained through the following relationship between the spontaneous emis- 
sion rates for an emitter in a homogeneous medium (FHom.) and into a waveguide (rwG)^ and 
the total classical dipole radiated powers in bulk (PHom.) and in a waveguide (/Vg) lfT3l : 

rwG ^ AVG ^^Q^ 
THom. Ph om. 

The validity of this expression is demonstrated in Appendix A (a similar relation for cavities has 
been established in 1241 ). This equivalence allowed us to use FDTD simulations to determine 
total spontaneous emission rates. 

3.3. Photoluminescence Collection 

With the help of Eq. (fTol i. we may describe the electric and magnetic fields radiated by a dipole 
located at z = inside the slab in terms of traveling supermodes. At an arbitrary z-plane along 
the waveguide direction, e = Lm=i ^/T,„/Tt„,e.x^{i^,nz) and h = ^ff^i ^r„,/rh„,exp(;^,„z), 
where e„, and h„, are the fields of the m-th eigensolution of Eq. with eigenvalue = 
/3„, + iou,m- The factor F,,, is the supermode emission rate from Eq. (|5]l. Applied in Eq. (|2]l, 
these fields lead to the following expression for the percentage of the total spontaneous emission 
carried by the fiber mode: 

rj,, = 2f /,„^+2 f W^A/^Re|/^exp[/(^,„-C)2]|, (H) 

m— 1 m..n=\ * ' k ) 



where 




IIsi^f^K)-zdS 



(12) 




'Re{//,(e,„ X h*J-zdS}Re[mef x h)) -zds} 
JJs{e,„xh})-zdS 



(13) 





and /„, = Re I \/f^^/J^ | • The factor of 2 is introduced to account for collection from both 

fiber ends. While the first sum in Eq. (fTTI) corresponds to the individual supermode con- 
tributions to r\pi, the latter corresponds to beating of the individual supermode field com- 
ponents. These contributions, which may be substantial, are periodic in z, with beat lengths 
Lz = 27r/(j3,„ — j8„), where j8„,,„ are the m-th and «-th supermode propagation constants. The 
m-th supermode individual contribution to the spontaneous emission collection efficiency is 
J?PL,m = fm • (r;n/r) = fmYm- The ratio = r,„/r givcs the fraction of spontaneous emission 
coupled into the supermode. Reflecting the fact that supermodes in our structure are hybrids of 
the fundamental fiber mode and slab TE and TM modes, /„, is called the fiber-mode fraction of 
the m-th supermode. 

The estimate of 77pL.,n could be improved through a more detailed consideration of the tran- 
sition between fiber-slab supermode and fundamental fiber mode. In the case of an abrupt 
transition (e.g., if the membrane is terminated at cleaved facets), modal reflection and trans- 
mission coefficients may be obtained through rigorous wave matching at the interface, yielding 
the power transferred to the isolated fiber mode ifTSl . On the other hand, a power transfer in 
excess of 90 % is estimated from the Fresnel reflection and transmission coefficients at such 
an interface, considering the supermode and fiber effective indices. Transmission levels of this 
magnitude were confirmed in FDTD simulations of a truncated slab of thickness 190nm, for 
dipoles placed at several distances from the transition. This suggests that the mode overlap 
method used here and in the previously presented FDTD results yields reasonable upper-bound 
estimates for the achievable collection efficiency. 

3.4. Supermode Calculation Method 

The eigenvalue equation Eq. (|4]i was solved with FEM with quadratic vector elements for Hj 
and Hy and Lagrange elements for H- (251. Given the symmetry of the problem, only half of 
the geometry was considered, with either anti-symmetric (x x E = 0) or symmetric (x x H = 0) 
boundary conditions on the y-axis, which automatically selected either hybrid-TE or hybrid-TM 
solutions. Perfectly-matched layers were defined over the outer boundaries to simulate an open 
domain, as leaky modes were expected. Perfect electric conductor (n x E=0) boundary condi- 
tions surround the PML on the outermost boundaries. The eigenvalue calculation produced a 
finite, discrete set of supermode solutions, which were used in Eq. (fTTT i to estimate tjpl. The 
respective emission rates were calculated with the guided-mode expression in Eq. (|5]). In 
all our calculations, fifty supermodes were used, which was sufficient to closely reproduce the 
FDTD results. 

3.5. Simulation results 

The PL collection efficiencies for jc-, y- and z-dipoles calculated with the supermode expansion 
method above are plotted in Fig. [5] The discrepancies in the y- and z dipole cases relative to the 
FDTD results are largely due to an insufficiently fine grid used in the FDTD simulations, where 
a full 3D simulation is run. In the FEM simulations, a cross-section of the structure is analyzed 
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Fig. 6. Electric field distribution for (a), (b): hybrid-TEQ types-I (slab-like) and II (fiber- 
like) supermodes; (d), (e): hybrid-TMo types-I and II supermodes. Fields are for f^iab = 
100 nm and A = 1.3 flm. (c) and (f): effective indices for (c) hybrid-TEg and (f) hybrid-TMo 
supermodes of types I and II for varying fsij,),- Fields have been normalized to the total field 
amplitude maxima in each case. Color scales are in arbitrary units. 



and hence can be discretized with a finer resolution. As discussed below, the main hybrid-TM 
modes, excited by either y- or z-dipoles, are strongly concentrated in the region between fiber 
and slab, requiring a very fine grid for high accuracy. To avoid prohibitively long computation 
times, a coarser grid was used. 

Hybrid TE modes: x-oriented dipoles 

Figure [7]; a) shows a representative effective index distribution for hybrid-TE supermodes. For 
A:-oriented dipole excitation, the two supermodes labeled I and II have the largest contribution to 
the total collected PL. Types-I and II supermodes are mostly concentrated in the slab and fiber 
regions respectively, and are thus referred to as slab-like and fiber-like. Representative electric- 
field profiles for these main supermodes are shown in Figs.|6ja) and |6jb), while Fig.|6jc) shows 
Re{«efT} as a function of fsiab for these modes. The periodic oscillaton of tjpl with z, observed in 
the FDTD simulations, is very well reproduced by the mode expansion calculations, as shown 
in the inset of Fig.|5ta). This oscillation is due to the beating of the two main supermodes, and 
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Fig. 7. (a) Effective indices for hybrid-TE supermodes at fsjaj, = 190nm and A = 1.3 flm. 
The two supermodes with highest modal emission rates are indicated with filled circles 
labeled I and II, corresponding to slab-like and fiber-like modes. Open circles indicate leaky 
modes, (b) Fiber fractions for supermodes I and II. (c) PL collection contributions from the 
individual supermodes, for an x-oriented dipole moment. The factor of 2 is to account for 
collection from both fiber ends. Black dots are for leaky modes. (d) Modal emission ratios 
for supermodes I and II. Filled squares/empty circles are for type I/II supermodes. 



persists when leaky eigenmodes are excluded from the calculation. The beat length is given by 
Lz = 2n/ (/3i — /3u), where jSiji are the respective propagation constants. 

The supermode power contributions, T7PL,m, are shown in Fig.|2lc). The slab-like supermode 
is dominant for fsiab ~ 160nm and fsiab ~ 240 nm, being surpassed by the fiber-like mode at 
intermediate thicknesses. In this range, the fiber-like supermode's significantly larger fiber 
fraction (Fig. |7jb)) is enough to compensate for its much lower modal emission ratio yo.ii 
(Fig- Eld)). Individually, leaky supermodes pointed out in Fig. 13 a) add little to the total col- 
lected power, however their aggregate influence is non-negligible: excluding leaky modes, the 
maximum collected power drops between 40 % and 70 %, a significant portion of the total 
collected PL. 

Among the leaky modes with TjpL.m above 1 % of the slab-like supermode contribution, 
the one with the lowest individual contribution has n^' S 156, corresponding to A^5=o.9 ~ 2.6 
wavelengths for 90 % amplitude decay. As a result, the highest efficiency collection shown in 
Fig. |5la) will be achieved for coupling lengths of a few microns; the contribution from leaky 
modes will begin to drop out as the coupling lengths increase beyond this. The slab- and fiber- 
like supermodes have A^g^g.g ~ 427 and A^g^o 9 ^ 22, respectively, so that a coupling length of ^ 
30/xm continues to allow collection of PL coupled to both main supermodes. Finally, collection 
that is at least as high as the slab-like supermode's contribution in Fig. |7ja) is achievable for 
lengths shorter than 550 jJ-m. We note that this latter result is still on par with what can be 
achieved through free-space collection (Fig.|5la)). 

Hybrid TM modes: y- and ;-oriented dipoles 

Figure [Sfa) shows a representative effective index distribution for hybrid-TM supermodes, 
taken at a slab thickness fsiab = 190 nm that is in the middle of the simulation range. For 
both y- and z- dipole excitation, two supermodes, labeled I and II, have the largest contribu- 
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Fig. 8. (a) Effective indices for hybrid-TM supermodes at fsjaj, = 190nm and A = 1.3 /im. 
The two supermodes with highest modal emission rates are indicated with filled circles 
labeled I and II. Open circles indicate leaky modes, (b) Fiber fractions for supermodes I 
and II. (c) and (e): PL collection contributions from individual supermodes, for y- and z- 
oriented dipole moments respectively (notice changing axis scales). The factor of 2 is to 
account for collection from both fiber ends, (d) and (f): Modal emission ratios, for y- and 
z-oriented dipole moments of type-I and type-II supermodes (notice changing axis scales). 
Filled squares/empty circles are for type I/II supermodes; black dots are for leaky modes. 



tion to the total collected PL. Representative electric-field profiles and Re{nefT} for these main 
supermodes are shown in Figs. lUd), |6te) and |6tf). The type-I supermode is well confined, 
having the lowest Imjnefr}, typically an order of magnitude lower than its type-II counterpart 
and all additional leaky modes. Leaky supermodes, pointed out in Fig. Uta), provide a very 
small contribution to the collected PL, typically at least two orders of magnitude lower than the 
supermode with the largest contribution to tjpl. The type-I supermode is mostly concentrated 
in the fiber for small fsiab, and in the slab at larger thicknesses, as apparent in its fiber frac- 
tion plotted in Fig. [8jb). The type-II supermode resides mostly in the fiber region for all slab 
thicknesses, although its fiber-fraction increases with fsiab- 

The individual supermode contributions tJpl.hi to the total PL collection for a y-dipole are 
plotted in Fig. [Sc). The type-I supermode contribution dominates at all thicknesses, while 
most leaky mode contributions range below 10^^. Its modal emission ratio is maximized for 
fsiab ~ 200 nm (Fig. [Sid)), for which a strong emission inhibition is observed (Fig. [3a)). In this 
regime, Yo,i is as high as 15 % of the total emission rate (30 % for emission in both ±z directions 
). The decreased type-I contribution at larger fsiab is due to a reduction of both its fiber-fraction 
/oj and modal emission rate 701, as shown in Figs.[8lb) and [Sd), respectively. Along with its 
increased fiber fraction /o n, this results in the type-II supermode power approaching that of the 
type-I supermode, in spite of its significantly reduced 70.11 ■ 

In comparison, the low T7pl obtained for z-polarized dipoles (Fig. [Sjc)) is a consequence 
of the small mode emission rates (Figs. [8je) and [Uf)). In both hybrid-TM and hybrid-TE 
modes, the z-electric field supermode components are dwarfed by the dominant, x- or y-field 
components, resulting in small effective areas for z-directed dipole moments. We note that the 
relatively poor coupling of z-dipoles can prove polarization sensitivity which can, for example. 



be used to determine the in-plane polarization of an embedded semiconductor quantum dot. 

As was the case for the jc-dipole, 77pl oscillates as a function of z, with a period that is dom- 
inated by the beating of these two main supermodes. In addition, tJpl decreases as a function 
of z, due to the decay of the various supermodes. However, this decay is significantly slower 
than that found in the x-dipole case, as leaky modes play a far less prominent role in PL col- 
lection here. Over all slab thicknesses, nj^ ^ 1440 for type-II supermodes, corresponding to 
Ns=Q.9 ~ 24 free-space wavelengths for a 10 % amplitude decay. A coupling length shorter 
than 31 jim thus allows for collection efficiencies at levels close to those plotted in Figs.|5](b) 
and HJc). For the type-I supermode, Ns=q.9 ~ 554, so that a collection level at least as high 
this supermode's individual contributions (Figs.Oc) and |Se)) may be achieved over coupling 
lengths of several hundreds of microns. 

4. Discussion 

The results of Section |3] indicate that laterally bound or quasi-bound hybrid supermodes play 
an essential role in luminescence collection in the single-emitter probing technique described 
in this work (and schematically summarized in Fig.[T]i. The overall PL collection efficiency de- 
pends on four interdependent parameters for each supermode: total (F) and modal (F,,,) spon- 
taneous radiation rates, fiber-mode fraction (/,„), and the supermode decay rate (oc^,,,,). High 
fiber-mode fractions may be achieved by phase-matching the fiber and slab modes, however, a 
power shift toward the fiber can result in a reduced field concentration at the dipole position, 
and therefore smaller modal emission rates. A compromise between F„, and /„, determines the 
optimal collection efficiency. The decay rate a-,,,, determines the length over which supermode 
collection is most efficient. 

A first method to obtain a substantially enhanced collection efficiency is through a smaller 
diameter fiber In Fig.|9ta), the maximum and minimum tjpl for an x-oriented dipole are plotted 
against slab thickness, for a 300 nm radius fiber (the same procedure which produced Fig. |5] 
was used, and the 500 nm radius results are re-plotted here). The maximum PL collection is in 
average 2.5 times that for a 500 nm fiber, while the minimum efficiencies are increased in av- 
erage by one order of magnitude, and vary between 2 % and 4 %. The collection enhancement 
is due to the increased slab-like and leaky supermode PL contributions (Fig.|9|b)), which result 
from increased fiber fractions (Fig.|9jc)) that are enough to compensate for a small decrease in 
the modal spontaneous emission rates. In addition, leaky modes play a more important role than 
in the R — 500 nm case, especially for fgiab ^ 180nm. Among the leaky supermodes with TjpL,m 
above 1 % of the slab-like supermode's, the one with the lowest individual contribution has 
nT^ ^ 167, or A^5=o.9 ~ 3 wavelengths for 10 % amplitude decay. As such, the maximum col- 
lection efficiencies in Fig.|9jb) will occur for coupling lengths that are within several microns. 
The slab-like supermode has A^5=o.9 ~ 382, so collection at least as high as its contribution in 
Fig Htb) is achievable over lengths shorter than « 500 /xm. 

Our calculations also show good prospects for probing InAs quantum dots in non-suspended 
slab configurations, where a AIq vGao.sAs (refractive index w 3.0 at A = 1.3 [im) substrate is 
present underneath the GaAs host layer. For a slab thickness of 250 nm, we find, similarly to the 
suspended membrane case, a slab-like supermode, several slab-concentrated leaky supermodes, 
and a fiber-like supermode, which in total yield a maximum collection efficiency of w 1.6 %. 
The reduced collection with respect to the GaAs membrane case is primarily due to a reduction 
in the supermode effective areas, as fields extend over longer lengths inside the substrate. In 
addition, the supermode decay lengths significantly change. The fiber-like supermode, which 
provides the largest contribution to PL collection, drops to 90 % of its original amplitude in 
approximately 2 /im. Probing lengths of a few microns are thus required for highly efficient PL 
collection. The slab-like and first slab-concentrated leaky supermodes have A^5=o.9 ~ 254 and 
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Ns=Q.9 ~ 24, respectively. For coupling lengths as long as w 3 1 ^m, the collection efficiency is 
at least 0.13 %, the contribution of these two supermodes alone. 

While the above cases and those studied in the rest of this paper largely involve a tradeoff 
between r,„ and /,„, situations exist in which spontaneous emission enhancement and high T7pl 
can be reached through small effective area, high fiber-mode fraction supermodes. For example, 
this possibility arises for the hybrid-TM modes if a y-oriented dipole is placed between the 
fiber and slab. At fsiab ~ 165nm, strong hybridization results in y-polarized supermodes with 
high confinement in the air gap, similar to the slot waveguide modes introduced in li26l . The 
spontaneous emission rate into this supermode alone reaches w 1 .7 times the homogeneous 
free-space rate, and the phase-matching between fiber and slab leads to a fiber-mode fraction 
close to 50 %. This indicates promising possibilities for fluorescence collection from dipoles 
located near the slab surface. Further investigation on single emitter spectroscopy appUcations 
utilizing this configuration is under way liZTl . 

For the case of jc-oriented dipoles, the collection efficiency may be effectively enhanced 
through stronger lateral confinement. In the hybrid-TE mode case of Section [3j PL collection 
efficiency is limited by the small fiber-mode fraction for the slab-like supermode and the even 
smaller fiber-mode fraction (and smaller spontaneous emission rates) for the slab-guided, leaky 
supermodes. To address these limitations, the slab may be converted into a suspended chan- 
nel waveguide, which has been fabricated in a number of material systems, including recent 
experiments involving semiconductor quantum dots li28l . This structure provides guided su- 
permodes with strong field confinement in the channel, so that the supermodes carry most of 
the emitted power, and may be designed to have high fiber-mode fractions, by phase-matching 
fiber and channel modes. As an example, for a channel thickness fchannei = 256 nm and width 
^channel = 250 nm, two bound modes exist and carry « 80 % of the total radiated power by 
an jc-oriented dipole at the channel center A total PL collection efficiency in excess of 70 % 
results. A detailed analysis of this system has been prepared li29l| . 

Finally, we consider the potential for resonant fluorescence measurements within our fiber- 
based probing scheme. As mentioned previously, resonant spectroscopy offers many advantages 
in comparison to non-resonant PL measurements, most notably improved spectral resolution, 
avoidance of generation of decoherence-inducing excess carriers and well-defined state prepa- 
ration m 121 [3] |4] |7] [12| • One of the main challenges in performing this type of spectroscopy 



is effectively isolating the residual or scattered resonant excitation beam from the resonance 
fluorescence signal, as they occur at identical wavelengths (such separation is easily accom- 
plished in non-resonant PL through a grating spectrometer or spectral filter). Recently, free- 
space optics collection of resonance fluorescence from a quantum dot inside a planar cavity 
was demonstrated |12|. In this work, separation of the fluorescence signal from the excitation 
beam was accomplished through an orthogonal excitation-collection scheme, where the sample 
was excited in-plane through a cleaved optical fiber, while the vertically emitted fluorescence 
was collected with a microscope objective. 

Through the setup of Fig. I2a), isolation of a resonant fluorescence signal can be achieved 
in a similar way. In a resonance fluorescence experiment, the pump signal (indicated in green 
in the figure) is tuned to the emitter transition. Illuminated by the resonant pump, the emitter 
radiates in both forward and backward directions of the fiber waveguide. The backward signal 
is routed to an external detector via a directional coupler, illustrated in Fig. 12 a), and is well- 
isolated from the pump signal (a fiber optic circulator could alternatively be used). Given the 
high fiber-collection efficiencies predicted in Fig. |5] (the collection in the backward channel 
will be O.Stjpl), this configuration may provide substantially superior signal-to-noise ratios in 
comparison to schemes that rely on free-space fluorescence collection. Alternately, we note that 
one might envision exciting the sample with a resonant free-space beam normal to the sample, 
while collecting the fluorescence through both the forward and backward channels of the fiber. 
This would create an orthogonal excitation-collection scheme similar to that of lfT2l . with an 
overall efficiency given by tjpl as plotted in Fig.|5] Such a configuration may be preferable in 
situations where reflections in the fiber path between source and emitter can not be avoided. 

5. Summary and Conclusions 

A technique based on the use of an optical fiber taper waveguide for collecting emission from 
single emitters embedded in thin dielectric membranes was analyzed with numerical simula- 
tions. The probing configuration was modeled as a composite waveguide formed by a micron- 
scale silica fiber on top of a dielectric slab. Finite-difference time-domain simulations were used 
to estimate the efficiency with which emission from in-plane and vertically-oriented dipoles can 
be collected into the fundamental mode of the fiber For in-plane dipoles, collection efficiencies 
superior to those achievable with a high numerical aperture objective are predicted, by as much 
as an order of magnitude. A maximum collection efficiency of ~ 7% may be obtained, with 
insignificant radiation suppression due to the slab. For vertically oriented dipoles, collection 
efficiency superior to free-space collection by at least an order of magnitude is predicted. A 
maximum collection above 35 % is obtained, albeit with strong slab-induced radiation suppres- 
sion (Purcell factor of < 0.05). It is important to note that the presence of the fiber does not 
significantly affect the radiative rate of the slab-embedded dipole, so the suppression equally 
affects free-space collection. 

Finite element simulations were used to understand the FDTD results through contributions 
from guided and leaky supermodes of the composite waveguide formed by the fiber and slab. 
Our analysis determined the essential role of laterally bound or quasi-bound hybrid supermodes 
in luminescence collection, as well as the relevant parameters for increased efficiency: total (F) 
and modal (F,,,) spontaneous radiation rates, fiber-mode fraction (/„,), and the supermode decay 
rate (Ot^m)- In short, the ratio r„,/r (or, the supermode j3-factor) and /,„ must be maximized, 
while a-,,n must be minimized for higher collection levels. Potential methods to further increase 
the collection efficiency, based on supermode analysis, have been introduced. Furthermore, 
probing of single quantum dots in unprocessed (non-undercut) samples is predicted to yield 
collection levels on the same order as obtained with free-space collection. Finally, we have 
described how the probing method is amenable to both non-resonant photoluminescence and 



resonant fluorescence measurements. 

In addition to efficient collection, optical fiber taper waveguides offer a number of benefits 
in the study of single solid-state emitters. These devices have standard single mode fiber input 
and output regions, and are therefore easily integrated with the technology of low-loss fiber 
optics. As discussed elsewhere Q, low-temperature measurement setups incorporating fiber 
taper waveguides and other fiber optic components have been used in cavity QED experiments 
in which precise knowledge of optical losses and power levels are a requirement. Furthermore, 
we note that as a movable probe, fiber taper waveguides have been used to interrogate two- 
dimensional arrays of devices on a chip II2TI . The results of this article indicate that these ad- 
vantageous features, previously demonstrated within the context of microcavity spectroscopy, 
can be utilized in direct single emitter spectroscopy. 



Appendix A - Classical point dipole radiation in the presence of a waveguide 

In this section, we demonstrate the validity of Eq. ( fTOl l. which states that the ratio of the spon- 
taneous emission rate in bulk to that in a waveguide is equal to the ratio of the total classi- 
cal dipole radiated power in bulk to that in a waveguide. We consider a point dipole source 
J(r,f) = — /op5(r)e^"**' radiating in the presence of an arbitrary, lossless dielectric waveg- 
uide extending along the z direction. The electromagnetic field may be expanded in terms of 
the bound and radiative modes of the waveguides as fOl : 



(A-1) 




The modes are such that the following orthogonality relations are valid: 
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The total power flowing through the waveguide at a position z > is obtained by integrating 
1/2 (E X H*) • z over the cross-section: 

|2 r r. (a\\2 ■ 
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We next normalize the guided and radiative mode fields as in Eqs. ^ and ([8]l. The guided mode 
propagation constants j3j may be expressed as lfT3l 
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Considering the waveguide material to be non-dispersive, the modal group velocity may be 
written as |[T3l 
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Using Eqs. (IA-7l i and iA-6i and Eq. ( IA-8b . the power through the waveguide cross-section, 
normalized to the homogeneous dipole power, Eq. ^ becomes: 



dp 
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The same expression is obtained by normalizing the spontaneous emission rate Y in Eq. dB-lOl l 
by the spontaneous emission rate of a dipole in a homogeneous dielectric medium, given by: 



Horn. 
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To verify Eq. ( fTOl i. we used FDTD to simulate a point dipole radiating at A = 1 .56 /zm near 
an optical fiber of index « = 1 .45 and radius a ~ 500 nm, surrounded by air The radiated power 
ratio Avg/^'hoiti. is compared to the spontaneous emission rate ratio Fwo/rHom., calculated 
with the semi-analytical expressions of iH. Figures [TOt a) and [TOl b) show the two ratios for 
radially and longitudinally oriented dipoles at position r along the radial direction. 

The relation between the classical guided fiber mode power P^/Py^a and spontaneous emis- 
sion rate F^,/Fwg ratios was also verified. The steady-state fields recorded at one of the waveg- 
uide edges were used in the overlap integral of Eq. (|2]i, together with the fundamental fiber 
mode fields, to yield Pg. Figure [TOf c) compares the ratios for the radially-oriented dipole as a 

function of its nnsitinn r in thf» radial Hirprtinn 

1.3^ : . ^ ^1 ^ ^ ^ 

(b) 
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Fig. 10. (a) Emitted power and modal spontaneous emission rate at A = 1.56^im, normal- 
ized to the corresponding quantities in a homogeneous medium, for a (a) radially and (b) 
longitudinally oriented dipole at position r in the fiber described in the text, (c) Emitted 
guided mode power, Pg, and spontaneous emission rate, Tg, normalized to the total emitted 
power and emission rates, at A = 1 .56 /Jm for a radially oriented dipole at position r in the 
fiber described in the text. 



Appendix B - Heisenberg equations 

To describe the spontaneous emission of a two-level emitter in the vicinity of a dielectric waveg- 
uide, we obtain the Heisenberg equations for the emitter interacting with a vacuum field reser- 



voir that is described in terms of the guided and radiative waveguide modes propagating in the 
Z-direction. 

Through the quantization procedure of 1301 , the positive-frequency electric field operators 



are defined as 
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Here, as in Section[3] fi = {f,in, co) labels a guided mode of order m traveUng in the direction 

/ •z, / = ±1, with propagation constant P^, and inverse group velocity = ^ = 
is the mode field profile and is the corresponding photon annihilation operator. In ( IB-21 ), 
V = {n,l5,(o) labels a radiative mode of order « and propagation constant P, with Cy the cor- 
responding field profile and Uy the photon annihilation operator The j3 -integral limits are such 
that evanescent radiation modes are excluded (e.g., |j3 1 < o/c if the waveguide is embedded in 
air) Ii23il . Also, Ny is such that 



dS (evxh;,)-z = A^v5(j3-j3') 
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The guided and radiative field profiles are normalized as in Section |3] such that [0^,0^/] = 
5{co - co')5ff5„,„,> and [ay, fly'] ^ 5{co - co')5{j5 -j5')5„„>. 

Next, consider a two-level atom (emitter) in the dipole approximation, interacting with a 
photon reservoir. The atom has a dipole moment deg for the transition with energy Iicoq between 
the ground \g) and excited \e) states. We define the atomic transition operators a = 
CT+ = |e) {g\, G. = \e){e\ — such that the atomic Hamiltonian is Ha = hcOoG-/2. The total 

Hamiltonian for the interacting atom and reservoir is 



H — Ha +Hr+Hi^. 

The interaction Hamiltonian for electric dipole coupUng is given by ///^ = - 
the rotating-wave approximation, may be written as [Sl: 

+'^nE^ t/fi) J dp 
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with 



{deg-t^) and Gy = 
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The formal solutions of the field operators are ||3T1 : 

2^,Ato) + Gly f dt' a-e'^"-'^^'' . 
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Under the Markoff approximations a (f — t) ^ O (t) and G^.v (o) « G^x^y (£00)5 the following 
Heisenberg equations for the atomic operators may be obtained: 



''^ -icoQO- - ^d- + a- \ Y,^^^/^gnUt) + Y, J V^ry,in{t)j , (B-8) 



dt 2 



^ = -r(i + (T,)-2a+|^^Vr;^^,,n(0+Eyf;rv,i„(f)j 

2a_ Vr;*4,„(r) +E / dp v^Vi.Jr)) , (B-9) 

\f.m n J J 

r= (Lr^+L/^J3rv), (b-io) 

\f,'n " •' ) 



with a = a-e''^\ 



^/r^=V2nG^{coo) and s/y\ = V2nGv{(Oo). (B-11) 
To arrive at these results, the following input field operators have been defined ||30l : 

'•v,in(0 J V27rJ-oo \ flv(fo) J 

Here, fl^.v(fo) = flM,v(fo)e'™o. 
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